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An automated linear measurement sys- 
tem (LMS) has been developed lo dc- 
t ermine the nonlinearity of a tuned 
30 MHz power detector over a 6.021 
dB range. This detector uses a single 
thermistor bead design with thermal 
isolation to obtain nearly linear tracking 
over a 4:1 change in input power. The 
nonlinear correction for this change, 
determined by the LMS, i:> on the or- 
der of 1.00030 (+ 130 jiB) for the de- 
tector pjcscntly in use. Initial 



experiments indicate an expanded tm- 
ccnainty of ±0.138% (±598 jiB), 
which is based upon TyjH; A and 
Type B components. 
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1, Introduction 

There has been a recent interest In and demand 
for a calibration service at NIST to support rf at- 
tenuators and voltage doublers that operate specif- 
ically at 30 MHz. The first step required to offer 
such a service is to develop a reference standard. 
For the best possible accuracy, a tuned single- 
element thermistor mount was chosen, A linear 
measurement system was designed and constructed 
at NIST to calibrate the nonlinearity of this mount. 
This paper contains a description of the LMS, an 
explanation of the measurement scheme, calibra- 
tion results, and an uncertainty analysis. 

2. System Description 

A diagram of the system is shown in Fig. 1. The 
signal generator provides a stable 30 MHz rf signal 
that is amplified and Filtered before the signal is 
split into two channels. One channel consists of a 
variable phase shifter and a fixed attenuator. A 
coaxial switch either terminates the signal with a 
50 fl termination or feeds the signal into a power 
divider, which splits the signal again. Half of the 



signal is detected by a single-element thermistor 
mount, Pi, and the other half is fed into the hybrid. 
The second channel consists of a variable attenua- 
tor. Like the first channel, a coaxial switch either 
terminates the signal with a 50 fl load or feeds the 
signal into a power divider, which splits the signal. 
Half of the signal is detected by a single-element 
thermistor mount, Pj, and the other half is fed into 
the hybrid. 

The hybrid takes the sum and difference of the 
two input signals. The difference is fed into a diode 
detector to rectify the signal, and then into a null 
meter. The sum is fed into a coaxial switch, and 
when the switch is in position 1, the signal is de- 
tected by a third thermistor mount, Pj, which is the 
thermistor to be calibrated. 

Each thermistor mount is connected to a NIST 
Type IV bridge and a digital voltmeter to measure 
rf power. 

The computer controls the signal generator, the 
digital voltmeters, and the switch controller, and 
handles the data acquisition and processing through 
an IEEE-488 bus. 
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Fig. 1. Block diagram of 30 MH2 linear measurement system. 



3. Design of 30 MHz 
Thermistor Mount 



Single-Element 



The dc response of detector Pj must be nearly 
linear with changes in input rf power. A single 50 fl 
thermistor bead design was selected as a linear de- 
tection scheme [1]. This detector is used in con- 
junction with a NIST Type IV self-balancing 
dc-substitution rf power meter modified to bias a 
50 n detector [2]. The T^pe IV power meter is de- 
signed to change the thermistor bead bias current 
so that the thermistor always maintains the same 
resistance. The detectors are placed in the LMS 
housing where the temperature is held to ±0.2 °C. 

It is difficult to filter the rf signal from the power 
meter leads due to the nature of the single bead. 
The total rf signal should appear across the bead in 
the ideal case. Other considerations in the detector 
design include thermal stability and any forms of rf 
leakage into or out of the detector. 

Several ideas are incorporated in the detector to 
eliminate these problems. Figure 2 shows the cir- 
cuit diagram of a single-element thermistor mount. 
LC filter sections are Inserted to filter the rf signal 



in the power meter leads. No ferrite material (such 
as a ferrite core inductor) is used near the thermis- 
tor, since ferrite components experience changes in 
impedance with changes in the rf power. These 
impedance changes lead to nonlinearities in the de- 
tector response. Therefore, air-cored inductors are 
inserted near the thermistor. Ferrite-cored conduc- 
tors are allowable in sections following the first fil- 
ter section because the rf power is sufficiently 
reduced and renders any impedance variation neg- 
ligible. The parallel LC filter is tuned to resonate 
at 30 MHz. A special thermistor-mounting struc- 
ture has been developed, and consists of an electri- 
cally insulated doughnut and two copper blocks, 
one on either side of the doughnut. The thermistor 
is placed inside the doughnut hole and is encapsu- 
lated in an air pocket by the addition of the copper 
blocks. The copper blocks also provide a large 
thermal mass so that the entire structure cannot 
experience rapid changes in temperature. The 
thermal time constant is much longer than the time 
required to perform a single measurement cycle. 
Double-sided copper-clad fiberglass boards are 
used in the exterior detector housing and in the 
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Fig. 2. Circuit diagram of singic-clcmcnf thermistor mount. 
AJl inductors are hand-wound and a!i capacitors are 1 (xF unless 
otherwise labeled. 

internal compartment walls which separate various 
filter sections. This type of construction and the 
use of capadtive feedthroughs in the compartment 
walls reduce the amount of rf leakage though the 
detector. These physical and electrical construction 
considerations provide a stable linear detector with 
low rf leakage. 

4. Measurement Methods 
4.1. Calibration of Pj 

Three power meters are used in theLMS— Pi, Pj, 
and P3. When a given power meter is read, the 
notation used is Pxvzry where X denotes the power 
meter (1, 2, or 3), F is 1 if channel 1 is switched on 
and if channel 1 is switched off, Z is 1 if channel 
2 is switched on and if channel 2 is switched off, 
and "r" denotes that this power is a reading and not 
a true value. Powers without a subscript "r" are true 
values. 

The calibration of Pj is achieved using a three- 
stage method. First, with switches 1, 2, and 3 each 
in position 1, enough power is applied so that a 
nominal 12 mW is incident on the three thermistor 
mounts. The phase shifter and variable attenuator 
are adjusted to balance the two channels, thus 
obtaining a null on the null meter. Readings are 
taken on Pi, P2, and P3 and are designated Piur, 
ftiir, and Psiir, respectively. Next, switch 2 is moved 
to position 2, so that power is only applied to the 
first channel, Reading are taken on Pj and Pj and 
are designated Pnoi and Pna,, respectively, A nomi- 
nal 12 mW will be incident on P| and approximately 



3 mW will be incident on P3. Finally, switch 2 is 
moved back to position 1 and switch 1 is moved to 
position 2, so that power is only applied to the sec- 
ond channel. Readings are taken on P^ and Pj and 
are designated Pjoir and Fjoir, respectively. A nomi- 
nal 12 mW will be incident on Pj and approximately 
3 mW will be incident on P3, 

The calibration constant of Pj, denoted C?)l, is 
calculated using 



Ch. 



"ill* 
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(1) 



The derivation of this formula can be found in 
Appendix A, The calibration constant is a measure 
of the detector's nonlinearity over the 6.021 dB 
power change, and is used as a multiplication factor 
to correct the ratio measured by P3, where 



Pirn 



— Lhl 



3]lr 



(2) 



4.2, Power Measurements 



The NIST Type IV power meter must be con- 
nected to an external dc voltmeter. The substituted 
dc power, Poi, is calculated from measured voltages 
using 



P<fc = 






(3) 



where K.ff is the output voltage with norf power ap- 
plied, Kun is the output voltage with rf applied, and 
Ru (SO il) is the resistance of the thermistor mount. 
Figure 3 shows the measurement sequence for a 
power calculation [3]. An initial K^f is taken; rf 
power is applied and [<,„ is measured; rf power ts 
removed and a final Votr is taken. The initial and fi- 
nal dc measurements are used with the V^^ mea- 
surement to calculate the power and correct for any 
mount drift, which is assumed to be linear. The cal- 
culated value of Knff in Eq. (3) is given by 



K« =K 



off.i 



h-ti 



(Yon.t" Kiff.i) J 



(4) 



where Kirf,i is the voltage reading taken before rf is 
applied at time /[, Vans 's the voltage taken after rf 
is removed at time ti, and r; is the time at which F™ 
is taken. 
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Fig. 3. Sequence for measuring power meter dc voltages. 

5. Results 

The calibration constant, Cm., which is a measure 
of the nonlinearity of mount Pj, was obtained by re- 
peated measurements. The average value of the 
three hundred trials taken so far is CEu. = 1.00030 or 
+ 130 ^.B. Long-term data are being accumulated 
to validate the calibration of P3. 

Table 1 shows a sample calibration. Powers read 
by the three thermistor mounts are displayed for 
each of the three stages. The calculated calibration 
constant is shown, along with the actual step in 
power for each leg of the system. 

Tabic 1. Sample calibration results of the 30 MHz Linear 
Measurement System's standard mount 



Switches I&2 


Switch 1 


Switdi 2 


P(I) 11.5194 


11.5089 


0.0000 


P(2) 11.4170 


0.0001 


11.3991 


P(3) 11.0907 


2.7644 


2.7756 


Leg 1; 6.0335 dB 






Leg 2: 6.0160 dB 






Ch : 1.00029 







6. Uncertainty Analysis 

6,1. Evaluation of Type A Standard Uncertainty 

Evaluationof a Type A standard uncertainty may 
be based on any valid statistical method for treating 
data. Examples are calculating the standard devia- 
tion of the mean of a series of independent observa- 
tions, using the method of least squares to fit a 
curve to data in order to estimate the parameters of 
the curve and their standard deviations, and carry- 
ing out an analysis of variance in order to identify 
and quantify random effects in certain kinds of 
measurements [4]. 



The calibration of the standard mount, P3, has 
been repeated three hundred times to determine 
the repeatability of the system. Tests were per- 
formed at various times of the day over several da^ 
to cover as many random factors as possible, includ- 
ing variations of environmental conditions and the 
operator's ability to renull the system. The sample 
standard deviation of the mean is ±0.000335% or 
±2 iiB. Long-term data are being accumulated to 
validate this calibration, and a control chart is being 
developed to track any possible outliers or drift. 

6.2. Evaluation of l^P^ B Standard Uncertainty 

Evaluation of a T^pe B standard uncertainty is 
based upon scientific judgment using all of the 
available relevant information. This includes previ- 
ous measurement data, manufacturers' specifica- 
tions, data provided in calibration reports, 
knowledge of the behavior of relevant instruments 
and materials, and uncertainties assigned to refer- 
ence data taken from handbooks [4], 

The Type B evaluation of standard uncertainty 
accounts for the following factors: 

a. Uncertainty in the dc voltmeter measurements. 

b. Uncertainfy in the Type IV power meters. 

c. Imperfect isolation between the two legs, 

d. Uncertainfy because Pxi^Pno, 

e. Effects of impedance changes in Pj. 

f. RF leakage. 

g. Spurious signals and harmonics, 

6.2.1 VoUmetcp Uncertainty The uncertainfy 
in the individual voltmeter readings may be deter- 
mined by taking the total differential of the power 
expression, Eq. (3), which gives 



dF=-^(K.(fdKoft-K« 



dK^) 



(5) 



The total differential of power, Eq. (5), may be de- 
termined by taking the differential of K*, Eq. (4), 
which gives 



dK«=(l-T)dK„ai + rdKrffj , 



where 



r= 



h-tj 
ti~ti 



(6) 



(7) 



The uncertainties, 6Voif,i, dK,tf.f, and dK«., in the 
measured values of VaH.i, K>ff,f, and Km, are based on 
the voltmeter manufacturer's specifications. 

Figure 4 shows the uncertainfy in the power 
measurement as a function of power level, assum- 
ing the powers are ratioed as they are in the Cfn, 
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calculation, Eq. (1). The power measurements, Fjnj 
and Piin, which are approximately 3 mW, result in 
uncertainties of 0.036%. The other power measure- 
ments, which are approximately 12 mW, result in 
uncertainties of 0.008%. The uncertaintj' of Cm. 
due to the voltmeter readings may be found hy in- 
serting the individual power uncertainties into Eq. 
(1), which gives 



Ck + ^v 



Pin-d,u * 



iJJ 



ftoi 


+ ^301 


■fjoi 


-^201 


P31I 


4-^311 



1/2 



Pr^x-Aiu 



(8) 



This results in an uncertainly of d[-= ±0.028% or 
± 122 nB. 



0.12 



l?l 




4 6 8 

Power (mW) 



12 



Rg, A. Power measurement uncertainty due to DVM when 
ratios arc taken. 

6J.2 Type IV Powcr-Mctcr Uncertainty The 
four fxjssible sources of uncertainties internal to the 
Type IV power meter are the reference resistors, 
the operational amplifier open-loop gain, input 
offset voltage, and input bias current. Larsen has 
shown that the uncertainties due to the Type IV 
power meters are negligible compared to those of 
the voltmeters [2]. 



6^3 Imperfect Isolation An uncertainty in 
the calculation of CfiLmay result if Pi and Pi are not 
zero when they are assumed to be. The thermistor 
mounts are square-law detectors and are not suffi- 
ciently sensitive to determine whether Pi and P2 are 
low enough in the "off condition to avoid signifi- 
cant errors. This uncertainty is derived in Appendix 
B. 

The corrected formula for CkLi taking imperfect 
isolation into account, is 



aL= 



PlUr 




|l + y| 



where 



and 



"lin 



yi = 



6b, 
bw\Qii 



(9) 

(10) 
(H) 



Here, yj is the measure of Pi's contribution to Ckl 
if it is not zero when it is assumed to be, and ^i is 
the measure of Pj's contribution to Cfji, if it is not 
zero when it is assumed to be. The terms, bna and 
fczui are the corresponding voltages to the powers 
Pi IIP and/'jii), respectively. Assuming/*) 10 aid /'toj are 
12 mW, hm and iaji are equal to 0.7746 V. The 
term, Qp,, is defined in Appendix B and its value is 
approximately -1. The ^'s represent the b\ 
which were assumed to be zero in Appendix A. If 
the isolation between the channels is 65 dB, as is 
stated in the power divider manufacturer's specifi- 
cations, then J6 1 = >1Z)2= 0.000436. Using nominal 
values for the P's and substituting the values into 
Eq. (9), an uncertainty of ± 0. 1 1 2% or ± 488 jiB is 
obtained, 

6.2,4 Uncertainty Because Pai'^Pii* An un- 
certainty in the calculation of Chl may result if 
p7*i\ '^/'.iin. The derivation of this uncertainty can be 
found in Appendix C 

The measure of nonlinearity, a, of thermistor 
mount three is 



1-0- 



»[, 



Chl Pu~P\u 



(12) 
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where Pu is the reading of mount three at the low 
level and Pi\i is the power reading of mount three 
at the high level. Assuming Chl= 1.0004 worst 
case, /'Lr=3 mW, and Pyb = \2 mW, a is calculated 
to be 3.334 x 10-^ The ratio Ky^/Km is given by 



K310 



1 + aPsw, 



If Psoir is 10% greater than Pnor, then 

Km = i^"(3xi-i«) =1.000010. 

^310 1 + Q:(3) 



(13) 



(14) 



The effect of this being nonunity results in a cor- 
rected formula for CjjL, where 



ClIL — 



PiiQi 



Pilir 



1/2 



^301 ^^31^ 
^31 Ir 



1/2 Tl 



L \ p. 



= 1.000005 



(15) 



using nominal \^lues for the P's. The result is an un- 
certainty of ±0.0005% or ±3 |xB. 

6.2.S Uncertainty Due to Impedance Changes in 
Pi An uncertainty in the calculation of Cm. may 
result if the impedance changes in P3 with power 
level. The derivation of this uncertainty can be 
found in Appendix D. 

The corrected formula for Chi., taking impedance 
changes in P3 into account, is 




Cm.— 



1~ Jlms TjiA (1 — I /310I ) 
1 — /lms Awl (1 ~ I -Tsiil ) 



(16) 



where Pm is the reflection coefficient of mount 
three at the high level, /sio is the reflection coeffi- 
cient of mount three at the low level, and /Ims is the 
reflection coefficient of the system looking into port 
three. 

The measured value of /310 was actually taken 
10 dB below Tjio instead of the 6 dB step, so the 
uncertainty should be conservative. Using nominal 
values for the P's and the measured reflection coef- 
ficients, the calculated uncertainty is ±0.014% or 
± 62 jj,B. 



Impedance changes in the diode detector with re- 
spect to power level are assumed to have a negligi- 
ble effect on the overall uncertainty. 

6.2.6 Uncertainty Due to RF Leakage It is dif- 
ficult to assign a quantitative uncertainty due to rf 
leakage. In order to reduce leakage, all coaxial ca- 
bles in the system were replaced with semirigid lines 
and, wherever possible, SMA connectors were used. 
For now, the uncertainty due to rf leakage is as- 
sumed to be ±10 M,B or ±0.0023%. 

6.2.7 Uncertainty Due to Spurious Signals and 
Harmonics After the signal is amplified and fil- 
tered, any harmonics are at least —92 dBc, while 
spurious signals are no greater than - 84 dBc. This 
will result in a measurement error that will affect 
the uncertainty of Chl. 

From Eq, (3), the calculated dc power is 



P^ = 20 (Vis -Vl^), 



(17) 



where Fdc is in mW. Assuming Kofr^l V, then 
F<,o = 0.63246 V for a calculated power of 12 mW 
and Von = 0.92195 V for a calculated power of 3 mW. 

A harmonic at —92 dBc results in an uncertainty 
of ±0.0025% for a voltage measurement. This un- 
certainty translates to power uncertainties of 
(12 ±0.0004) mW and (3 ±0.0008) mW. Using the 
nominal powers to calculate Chl, an uncertainty of 
±0.023% or ±100 ]lB is obtained. 

A spurious signal at — 112 dBc results in an un- 
certainty of ±0.00025% for a voltage measurement. 
This uncertainty translates to power uncertainties 
of (12 ±0.00004) mW and (3 ±0.000085) mW. Us- 
ing the nominal powers to calculate Cm., an uncer- 
tainty of ±0.0025% or ±11 m,B is obtained. 

6.2.8. Overall l^pe B Uncertainty For each 
Type B component, an estimated range, ±aj, is 
given, assuming that the quantity in question has a 
100% probabihty of lying within that interval. The 
quantity is treated as if it is equally probable for its 
value to lie anywhere within the interval. Therefore, 
it is modeled by a rectangular probability distribu- 
tion. The best estimate of the standard deviation, 
Uj, is 



M;=— -^ 



V3 



(18) 



Table 2 show^ all of the Type B components along 
with their corresponding uncertainties and standard 
deviations. The overall standard deviation of the 
Type B components, calculated using the root-sum- 
of-squares method (RSS), is ±0.069% or ±299 m,B. 
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Table 2, Type B components of the 30 MHz Linear Measure- 
ment System with corresponding uncertainty ranges and stan- 
dard deviations 



Component 


Range 


Standard deviation 




(7c) 


i%) 


dc voltage measurements 


±0,028 


0.016 


Imperfect isolation 


2:0.112 


0.065 


Pxi^Pjio 


±0.0005 


0.0003 


Impedance changes in Pj 


±0.014 


o.ooa 


rf leakage 


i0.0023 


0.0013 


Spurious signals 


s 0.0025 


0.UO14 


Harmonics 


±0.0B 


0.013 


Combined Type B Standard 






Uncertainty (RSS) 




± 0.069 

( + 299 jiB) 



63 Combined Standard Unccrtaint>' 

The combined standard uncertainty, u^, is taken 
to represent the estimated standard de\iation of the 
result. It is obtained by combining the individual 
standard deviations, u-,, whether arising from a Type 
A or a Type B evaluation [-1 j. The technique used to 
combine the standard deviations is the RSS 
method. 

The total uncertainty reported is the expanded 
uncertainty, U, which is obtained by multiplying u^ 
by a coverage factor, /c. To be consistent with cur- 
rent international practice, the value of k used at 
NIST for calculating Ui$k~2. The total expanded 
uncertainty, U, of the mount's nonlinearity is calcu- 
lated to be ±0,138% or ±598 \lB. 

7. Conclusion 

The first step toward offering a calibration ser- 
vice at NIST to support rf attenuators and voltage 
doublers that operate at 30 MHz has been com- 
pleted. The reference standard, a tuned single- 
element thermistor mount, has been caHbrated us- 
ing an LMS, designed and constructed at NIST. 
The next step is to modify the LMS so that a device 
under test may be inserted into the system and cal- 
ibrated against the reference standard. 

8. Appendix A. Scattering Coefficient 
Analysis of the LMS 

The 30 MHz Linear Measurement Sj"stem can be 
considered a five port network, as shown in Fig. 5. 
The analysis assumes that the LMS is an ideal sys- 
tem—that is, the system is linear, and when hi and 
6; are assumed to be 0, they truly are. 

The pertinent scattering coefficients are 



l>i=S2it2i + 5^:^12+ Sua i + Si^o 4 + S2iiii 
bs — Sjia 1 +5i!0 2+ 551^3 + Ss4fl4 + Si5fl 5 . 



(19) 



If(it = /'i6i and<ii=n6j, then Eq. (19) may be writ- 
ten as 

b^{l-r\Sn) = Sur2bi-k-Suai+S»ax+SKai 

bli 1 - flSn} = Sjl / "1 fcl 4- 523 as + ^24 04 + ^25 os 

bi=Sii Fibi + Sji /2&1 + 5iia3 + Sm04 +83505 
bi = Stti'ibi +St2ribi + S*^aj+S4iat+S4sai 
b5 = S§i I'lbi + Ss2r'2bz + Siiai + Si404+ Sssos 

(20) 

This gives five equations and eight unknown vari- 
ables— aj, a*, flj, bu bz, bit bi, and ftj. Four variables 
can be eliminated from Eq, (20), say at, bt, cj, and 
bi. Eq. (20) becomes 



b, = P\:a3 + Q'nbz + Qhbi, 



(21) 



where the primed quantities are functions of the 
scattering coefficients and Fi and A. 

It is assumed that F^ is constant with respect to 
varying power levels, If now a} = Fjbit 

bi = iPhF3 + Qi3)bi + Ql2bz^ 

Qiibi + Qnbi (22) 



or 



, _ bi-Qnbi 
0.1 - — 1~ 

^13 



(23) 



where two complex constants, Qn and Q 13, describe 
the relationship among bi, b^, and bi. 
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Rg. S, LMS port nomenclature. 
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